Entangling the lattice clock: Towards Heisenberg-limited timekeeping 
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We present a scheme for entangling the atoms of an optical lattice to reduce the quantum pro- 
jection noise of a clock measurement. The divalent clock atoms are held in a lattice at a "magic" 
wavelength that does not perturb the clock frequency - to maintain clock accuracy - while an open- 
shell J = 1/2 "head" atom is coherently transported between lattice sites via the lattice polarization. 
This polarization-dependent "Archimedes' screw" transport at magic wavelength takes advantage 
of the vanishing vector polarizability of the scalar, J = 0, clock states of bosonic isotopes of divalent 
atoms. The on-site interactions between the clock atoms and the head atom are used to engineer 
entanglement and for clock readout. 
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Quantum entanglement is a crucial resource in quan- 
tum computing and has the potential to improve preci- 
sion measurements i^] . Here we propose a scheme for 
entangling an ensemble of several thousands of neutral 
atoms, with the specific goal of demonstrating the power 
of entanglement for measuring time. 

Measuring time with atoms relies on the fact that the 
quantum-mechanical probability of making a transition 
between two clock levels depends on the detuning Aw of 
the probe field lo from the atomic transition frequency wq- 
By measuring the probability as a function of uj, one can 
infer if the two frequencies are equal and thereby "lock" 
a local oscillator to the atomic transition. Counting the 
number of oscillations of the local oscillator tells time. 
The precision of measuring Aw is limited by the quan- 
tum projection noise 0]. For a measurement of N unen- 
tangled atoms the resulting signal-to-noise of Aw scales 
as \/]V: the standard quantum limit (SQL). The use of 
entanglement holds the promise of improving clock preci- 
sion to the Heisenberg limit, with signal-to- noise scaling 
as N. 

Measurements with uncertainty below the SQL may be 
achieved with squeezed atomic states |4j, 5, 6]. While this 
technique can address large number samples of atoms, 
squeezing experiments have attained signal-to-noise ra- 
tios far from the Heisenberg limit. In other work, clock 
measurements at the Heisenberg limit have been demon- 
strated for small numbers of entangled ions in traps 
0. Those experiments created maximally entangled 
states to achieve measurements at the Heisenberg limit 
via a Ramsey-type measurement protocol: a generalized 
7r/2 pulse creates a Greenberger-Horne-Zeilinger (GHZ) 
state, the atoms undergo free evolution, and a final gen- 
eralized 7r/2 pulse is used for readout. However, scaling 
ion traps up to entangle increasingly larger number of 
ions currently remains a work in progress. 

Using atom-atom interactions to engineer entangle- 
ment between neutral atoms trapped in a lattice may 
offer the best of both worlds: maximally entangled 
large-number samples. Previous proposals have noted 
the virtues of using alkaline-earth-like atoms in lat- 
tices for quantum information and quantum computing 



d, [13, [O, [13 ; here we focus on entangling atoms for 
improving the atomic clock. 
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FIG. 1: (Color online) Schematic of the entanglement pro- 
cess, (a) With (f> = Q a. single head atom (orange circle) and 
several clock atoms (blue circles) are trapped in the minima 
of a 1-D optical lattice, with one or fewer atoms per site. Due 
to an intensity differential of the underlying lattices, the clock 
atoms couple strongly to the cr+ lattice (solid blue line). The 
head atom is placed in a superposition of atomic states: one 
which couples strongly to the <j+ lattice and one which cou- 
ples strongly to the (t_ lattice (dashed orange line), (b) As (j> 
is increased, the latter state spatially separates and is trans- 
ported along the lattice, (c) This portion of the head atom is 
then brought into contact with a clock atom to entangle the 
two atoms, (d) The head atom is transported further to ob- 
tain entanglement with the remaining clock atoms in a similar 
manner. 



In optical lattice clocks, millions of divalent atoms 
(such as Sr or Yb) are trapped in an optical lattice (regis- 
ter) operating at a "magic" wavelength Am . At this wave- 
length, both clock levels, ^5*0 and ■^Pq, are shifted by the 
lattice lasers equally, so that the clock frequency remains 
unperturbed [13]. These clocks have demonstrated long 
coherence times and have already realized great improve- 
ments in both accuracy and precision over the current 
primary frequency standard [l4 . [l5| . Notably the signal- 
to-noise ratio of current-generation clocks is approaching 
the SQL [31, so an "entangled clockwork" may be of 
practical benefit. 
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We would like to entangle a string of strontium clock 
atoms. Each Sr atom occupies an individual lattice site 
as shown in Fig. [TJ To maintain the clock accuracy 
we require all lattice fields be at the magic wavelength 
of strontium. Entangling the atoms using short-range 
atom-atom interactions requires the transport of atoms 
between lattice sites. While coherent transport in opti- 
cal lattices has been demonstrated before using lattice 
polarization p^ . [l9| . such techniques would not 

work for the clock states of bosonic isotopes, which have 
a scalar (F = J = / = 0) nature: the optical potential 
does not depend on the polarization (be it circular or lin- 
ear) of the lattice lasers. Instead, we use a single J — 1/2 
"head" atom, which is transported from site to site using 
the optical polarization f23|. 

Lattice — Consider the superposition of two spatially- 
displaced standing waves of opposite circular polarization 
(d-±). The resulting E- field reads 

27r 27r 
E(z) = E+cr+cos{- — z) + £'_(t_cos(-- — z — ip), (1) 

where winding the phase similarly to the "Archimedes' 
screw" controls a relative displacement between the nodes 
of the two standing waves. The resulting optical potential 
reads 
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U{z) = C/o+cos^(— z) + U^cos\—z - ip). (2) 
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For an atom in an \F,Mp) state (with the quantization 
axis taken to align with the lattice lasers) 



^\ UM±^aUu^,^)]. (3) 



Here Q^{LiJm) and ap(wm) are frequency-dependent scalar 
and vector (axial) polarizabilities. Neglected tensor con- 
tribution is suppressed [21] for J — 1/2 atoms. 

The two clock states |0) = \^Sq) and |1) = pPo) 
will experience the same trapping potential; at A^ the 
two ac polarizabilities are the same. Note that the vec- 
tor part of the polarizability is zero for the scalar clock 
states. If two states of the head atom | t) = |F, Mp) and 
I I) — \F' ,M'p) have different vector polarizabilities, they 
will see different potentials. For an appropriate choice of 
lattice parameters, discussed below, the | [) state couples 
preferentially to the cr"*" permanent sub-lattice, while | ]) 
couples preferentially to the moving a~ sub-lattice. Un- 
fortunately, the commonly employed lin Z lin transport 
lattice algorithms [l^, [ij, [l^, cannot be employed 
here directly, as the potential will wash out for the clock 
atom in the lin _L lin configuration. However, for appro- 
priate choice of lattice intensities (discussed below), the 
clock atoms remain pinned to the (T+ permanent sub- 
lattice. This state-selective transport enables the entan- 
glement of the clock atoms with the head atom state, as 
shown in Fig. [T] and described below. 

Clock protocol — We describe the system as the prod- 
uct of the state of the clock register and the state of 



the head atom. For example, a possible basis state 
of TV = 3 clock atoms and a head atom in the spin- 
up state is = |110)| t) = |6)| ])■ We require two 
gates: a single-qubit Hadamard gate H (an analog of 
a 7r/2 pulse) and a two-qubit phase gate. No indi- 
vidual addressing is required. The phase gate Pi in- 
volves the "head" atom state-selectively transported to 
overlap with the target clock atom at position i in the 
lattice register, P^]..., 0^, ...)| t> = 0^, ...)| t) and 

p,|...,i„...)| T) = -|...,i.,...)| T). 

A practical realization of the Hadamard gate involves 
interaction with a near-resonant pulse of optical fre- 
quency oj for the clock qubits and a near-resonant pulse 
of microwave frequency lo' for the head atom. We want 
to measure the clock frequency ujq by tuning the driving 
frequency uo. Below we show that, as in the conventional 
Ramsey-type clock frequency measurement, the probabil- 
ity of making a clock transition depends on the detuning 
Aw = 07 — 070, allowing to zero-in on the clock frequency. 
Theoretical analysis is simplified by transforming into a 
rotating reference frame; the relevant chain of operators 
for the compound wavefunction involves a product of N 
clock-state rotation operators at frequency uj and a rota- 
tion operator for the head atom at frequency 07'. 

We start by filling a 1-D lattice with a single head atom 
and N clock atoms, = |000...)| [). Next we apply the 
Hadamard gates to the head atom and to all the clock 
atoms: 
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where \p) ^ |0...00), |0...01), is the computa- 
tional (binary) basis set for the clock register. In the 
next step we use the transport lattice and move the | t) 
state of the head atom along the clock register to perform 
a collisional phase gate at each site. 

We presume that the |0) and |1) clock states will have 
different scattering lengths for their interaction with the 
I ]) state of the head atom, and thus a different mean- 
field interaction. The transport lattice will be moved 
in such a way that the head atom will remain on-site 
with each clock atom for a sufficieiit period of time to 
produce a relative phase shift of tt This generates 

the entangled wavefunction 



1*2 



Here fcp — pj, p — X^jLo^ -Pi 2-', i.e., it is the number of 
raised bits in the binary representation of p. By applying 
the Hadamard gate to the clock register we obtain the 
GHZ state 
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(|000... 000)1 i> + |lll... 111)1 t)) 



The entire procedure may be considered as a generalized 
7r/2 pulse in the space of maximally-polarized states. 
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As in traditional Ramsey spectroscopy we let the GHZ 
state evolve freely for a time T . In the rotating reference 
frame, a phase x — (NAuj + Alu') T is accumulated dur- 
ing the free evolution. Notice that Aco is the detuning 
for the optical clock frequency, while Aw' is the detuning 
for the microwave transition of the head atom. There- 
fore Auj' ^ NAlli. Finally, we repeat the generalized 7r/2 
pulse, arriving at 

|*fi„ai) = |0...0){cos| II) -J sin I IT)}, 

where the A^-enhanced phase is encoded into the state of 
the microwave head qubit. 

For readout of the clock, the state of the head qubit 
can be read directly via laser-induced fluorescence. In 
the event that the light collection efficiency is insufficient 
to read out the state with high efficiency ^] , the state of 
the head atom can be transferred to many alkaline earth 
atoms (as done in the first stage of the entanglement 
algorithm) for efficient readout. 

Choice of the "head" atom and lattice parameters — 
We start with a detailed analysis of the transport lattice 
which dictates the choice of the head atom. We introduce 
the fractional intensity misbalance for the two circularly 
polarized sublattices of Eq. ([T]): S = {El - E'^)/{E\ + 
E"^). It does not depend on atomic properties and \5\ < 
1. For the clock atoms, the minimum depth of the optical 
potential — which occurs at (/? = (n-f l/2)7r, or positions 
(b) and (d) in Fig. [T] — is proportional to \S\. However, 
the larger \5\ is, the weaker the lattice becomes for one of 
the head atom states. The fixed value of (5 = -1-1/4 will 
be used in all the following calculations. 

We further introduce a ratio of the vector and scalar 
polarizabilities, 

The ability to translate the |t) state while holding the |J,) 
state stationary is determined by this ratio. For example, 
for p| — —1, Uq = and the atom does not see the 
stationary (t+ lattice. Similarly, for p| = +1 the atom 
has no coupling to the moving a~ lattice. In general, for 
a positive value of S we must satisfy the criteria — 1/(5 < 
Pi < —S and pi < 1/S together with pi > 5. 

Qualitatively, for the head atom transport one needs 
|a|n(ci;m)| ~ |a^(ti-'m)|. It turns out that at the rnagic 
wavelengths specific to divalent atoms (see Ref. [23|). 
none of the commonly used alkali-metal atoms satisfies 
this constraint. Alkali atoms have the nS 1/2 ground 
states and the vector polariz ability (rank 1 tensor) arises 
only due to relativistic effects. Fortunately, atoms with 
nPi/2 ground states (group HI) have large vector polar- 
izabilities. Aluminum is a suitable choice for the head 
atom. 

^^Al has a 3pi/2 ground state. The nuclear spin of 5/2 
gives rise to two hyperfine structure levels: F = 3 and 
F = 2, separated by 1.5 GHz. Cooling Al has already 



been demonstrated [24| with the goal of atomic nanofab- 
rication. The laser cooling was carried out on the closed 
3p3/2 ~ 3(^5/2 transition with the recoil limit of 7.5 pK. 
Once trapped, the atoms can be readily transferred from 
the metastable 3p3/2 cooling state to the ground (head) 
state. Lattice-trapped Al was also considered for quan- 
tum information processing j25||and for a microwave lat- 
tice clock (niicroMagic clock) |2l| . 

To evaluate the dynamic polarizabilities for Al, we em- 
ployed ah initio methods of relativistic many-body the- 
ory. To improve on the positions of atomic resonances, 
for low-lying energy levels we replaced the ab initio en- 
ergies with experimental values. The resulting dynamic 
polarizabilities of Al are shown in Fig. [21 
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FIG. 2: (Color online) Dynamic polarizability of Al as a 
function of lattice laser frequency. All values are given in 
atomic units. Marked points on the plot correspond to magic 
wavelengths for clock transitions in divalent atoms 23, 26]. 

We focus on the promising case of the Sr "blue" 
magic wavelength [26l| at = 389.9 nm. Here the 
clock atoms are confined to minima of laser intensity, 
reducing photon scattering. Scalar polarizabilities are 
asr(^m) — -470 a.u. and a\^{ujm) — -340 a.u. ^"^Al of- 
fers many viable |F, Mp) states to implement transport: 
IT) = |2, 2), |3, -3), or |3, -2) {p « -0.84,-1.25, -0.84, 
respectively) and |i) = |2, -2), |3, 3), or |3,2) {p w 
0.84,1.25,0.84), for example. We choose |t) = |3, -3) 
and li) = |2,-2). 

Due to being "blue-detuned" for both Al and Sr, 
the atoms will be confined to regions of intensity minima. 
While the transport lattice provides axial confinement, 
the radial confinement is provided by two transverse 
"blue-detuned" magic lattices with polarization parallel 
to the transport lattice's k-vector (to avoid interference 
with the transport lattice). The transverse blue lattices 
create a series of tubes [27|, each containing a transport 
lattice. This provides tight radial confinement and en- 
ables many transport lattices to be run in parallel. We 
note that to prepare the initial states in the resulting 3-D 
lattice, both species can be prepared in the ground vibra- 
tional states of individual lattice sites using 3-D Raman 
sideband coohng 
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The depth /S.U of the optical potentials is critical, 
as it determines the rate of unwelcome diffusion out of 
sites [1]. The depth varies with the displacement phase 
(f. For our choice of | t)j the depth is weakest when it 
is on-site with Sr. At this maximum-overlap position, 
A;7(|T)) = 2lL^/c\al,{Lo„,)\{l+5p^) and A[/(clock) - 
2/L7r/c|a|j,(wm)|, where II = f;; {E\ + E'i) and c is 
the speed of light. We require AJ7 > bEn, where 
Er — {2nh/X,n)'^/{2M) is the recoil energy for an atom 
of mass M. This translates into a minimum intensity of 
II ~ 20kW/cm^, determined by the lighter Al. 

Gate times and decoherence — The number of clock 
atoms which may be entangled will be limited by deco- 
herence and phase-gate times. The main sources of deco- 
herence are anticipated to be inelastic collisions between 
the clock atoms and the head atom and light scattering. 

The decoherence rate due to the photon scat tering is 
r-i = rj§^u;l \a%cj^)f II- Here « l/2^En/AU 
is a suppression factor [8,] accounting for atomic wave 
functions being centered at zero intensity. We find ~ 
10 s for Sr and « 8 s for Al. 

To estimate the time required for a phase gate op- 
eration, we note the interaction energy of two par- 
ticles in overlapped ground states of independent 3- 
D anisotropic harmonic potentials is given by 6E = 

Ili=xyz i™^)]''^ where ascatt is the scattering 
length, TO is the reduced mass, and (toZJ)^ is defined anal- 
ogously to the reduced mass with u)i being trap frequen- 
cies. 

The axial trap frequencies may be determined from the 
parameters above (/l, (5, etc.) at overlap. The radial trap 
frequencies are determined assuming the transverse lat- 
tices are operated at the same intensity as the transport 
lattice. Estimating the difference in the excited-state and 
ground-state scattering lengths to be ascatt ~ 100 a.u. 
yields an estimated gate time of r 20 /is and a trans- 
port time of T « 10 /iS. 

Conclusions — With the pessimistic assumption that a 
single photon scattering event will decohere all the clock 
atoms, we expect that - within each 1-D lattice - one 



would be able to put ^ IQ^ Sr atoms into the maxi- 
mally entangled GHZ state with high probability. This 
would enable a reduction in the projection noise of lat- 
tice clocks. Moreover, we note that many of the usual 
requirements for producing highly-entangled states be- 
tween atoms - such as single-site addressability, single- 
site readout, and unity site occupation - are absent in 
this scheme. 

Possible improvements would include the use of entan- 
gled states less sensitive to photon scattering as a source 
of decoherence. Moreover, it is possible that combining 
these techniques with more sophisticated gate operations 
could lead to the development of a full quantum com- 
puter. In that case, error correction techniques could po- 
tentially be used to further increase the number of atoms 
while maintaining high fidelity. 

But even without these improvements, we note that 
this scheme occupies an interesting "middle-ground" of 
experimental schemes for clock entanglement. It holds 
promise for use with larger numbers of atoms than has 
been demonstrated to date with ion traps. And while it 
cannot entangle as large-number samples as are used in 
spin-squeezing experiments, it may be able to produce 
greater levels of entanglement. 

Unanswered questions remain, such as the value of 
the Al-Sr scattering length (which determines the gate 
time) and the rate coefficient for inelastic Al-Sr collisions 
(which is an additional source of decoherence). With 3 
naturally occurring bosonic isotopes of Sr, it is likely that 
a favorable scattering length can be found. The inelas- 
tic rate is unknown, but we note that measurements of 
Al-group atoms in the ground state have observed 

slow hyperfine relaxation in collisions with J = noble 
gas atoms [soj . 
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